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The CoUatz, or 3a; + 1, Conjecture claims that for every positive 
integer n, there exists some k such that T''{n) = 1, where T is 
the CoUatz map. We present three cellular automata (CA) that 
transform the global problem of mimicking the CoUatz map in 
bases 2, 3, and 4 into a local one of transforming the digits of 
iterates. The CAs streamline computation first by bypassing cal- 
culation of certain parts of trajectories: the binary CA bypasses 
division by two altogether. In addition, they allow for multiple 
trajectories to be calculated simultaneously, representing both 
a significant improvement upon existing sequential methods of 
computing the CoUatz map and a demonstration of the efficacy 
of using a massively parallel approach with cellular automata to 
tackle iterative problems like the CoUatz Conjecture. 

Keyphrases: CoUatz Conjecture, massively parallel, determin- 
istic computational model, cellular automata 

1 INTRODUCTION 

Consider the CoUatz map T : M M defined by 




mod 2. 



mod 2 
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Given n, define the total stopping time to be tiie smallest k for which T'^ (n) = 
1. The Collatz Conjecture claims that every n has a finite total stopping time; 
it has been verified for inputs as high as 20 • 2^^ by Oliveira e Silva [IJ, 
and a generalization of the Conjecture by Kurtz and Simon was shown to be 
recursively undecidable [2]. If we define the stopping time to be the smallest 
k for which T^{n) < n, the claim that every n has a finite stopping time is 
equivalent to the Collatz Conjecture. 

The trajectory for an input n is the sequence 

n,T{n),T^{n),T^{n),... 

We can categorize the possible behaviors of the trajectory as follows [3]: 

(i) Convergent trajectory: T*^ (n) = 1 for some positive integer k. 

(ii) Non-trivial cyclic trajectory: T^(n) becomes periodic and T^{n) ^ 1 
for any positive k. 

(iii) Divergent trajectory: lim T^(n) = oo. 

The iterative nature of the problem and the complexity of behavior that 
arises from the simple premises of the Collatz Conjecture suggest the use of 
cellular automata to mimic computation. 

Notable instances of using deterministic computational models to simu- 
late the behavior of the Collatz map include De Mol's 2-tag system using a 
set of three production rules over an alphabet of three symbols to mimic the 
Collatz map [4], Michel's two one-tape Turing machines [5] whose halting 
problems depend on generalizations of the Collatz problem, and Bruschi's 
two one-dimensional cellular automata that calculate trajectories using tag- 
ging operations to detect parity [6]. In this paper, we demonstrate one three- 
dimensional and two two-dimensional cellular automata that also mimic the 
behavior of the Collatz map but whose evolution laws operate independent of 
parity, depending solely on local relationships between digits of the iterates. 
We find that our automata have the potential to compute trajectories more 
efficiently than do current methods both by bypassing computation of large 
parts of trajectories and by verifying an arbitrary number of inputs in parallel. 

2 CAl- BASE 3 

This CA acts upon an three-dimensional grid of cubical cells in two layers in- 
dexed by and 1, respectively. Each cell in the bottom layer is in one of four 
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states represented in Figure |2] by the corresponding colors in parentheses- 

(light gray), 1 (gray), dark gray, and empty (white) — and each cell in 
the top layer is in one of four states also represented in Figure |2] by the cor- 
responding colors in parentheses- odd-normal (gray), odd-special (odd- 
normal with a cross), even (light gray), and unknown-parity (black). In the 
planes of each layer, rows are indexed top to bottom by IN U {0} and columns 
right to left by Z, where row i and column j of the bottom layer lie respec- 
tively beneath row i and column j of the top layer. We denote a cell in row 

1 and column j of layer k by the ordered triple {i, j, k) and its state at time 
period t by Sij-,fc(t). 



As shown in Figure la given a cell in the bottom row (dark gray), the 



neighborhood of that cell is the set of cells labeled B-F (light gray), where cell 
B lies directly under cell F. As shown in Figure [lb] given a cell in the bottom 
row (dark gray), the neighborhood of that cell is the set of cells labeled D, G, 
H (dark gray). In other words, the neighborhood of cell 0) is the set of 

cells{(z-l,j,0),(z-l,j- l),(^-l,j-l,0),(^-l,J-l,l),(^,J-l,0)}, 
while the neighborhood of cell 1) is the set of cells {{i,j,0), — 

1,1)}- 




(a) Bottom layer (b) Top layer 

FIGURE 1 : Neighborhoods for cells in different layers 



If input n has ternary representation J2i=Q 3*^1, the CA is initialized as 
follows: for i from to iV and an arbitrary k, we have so,fc+j,o(0) — di. All 
cells in the top layer are set to state unknown-parity, and all other cells in 
the bottom layer are set to state empty. 

For simplicity, the rules in Figure|2]are presented in two dimensions: rules 
of categories la through Ic represent transformation of the bottom layer, 
while rules of categories Id through Ig represent transformation of the top 
layer. 
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FIGURE 2: CAl Evolutions Laws 



A rule in category la, lb, or Ic is represented as a 2 x 2 grid with an arrow 
directed at a 2 x 1 grid (for convenience, the 2x2 grids of rules in Figure|2]that 
involve the same 2x1 grid are grouped together). They correspond, respec- 
tively, to the bottom and top halves of the configuration in Figure [Ta] Specif- 



ically, in the 2x2 grid, the bottom-left square corresponds to A in Figure la 



the cell that is evolving, and its color denotes the future state into which that 
cell evolves given the states in its neighborhood. If this cell has coordinates 
(z, j, 0), proceeding clockwise, the grid's other cells correspond to the current 
states of cells B{i — 0), C{i — 1, 0), and D(i, j — 1, 0) respectively. 
In the 2x1 grid, the left and right cells correspond to cells F{i — 
and E{i — — 1, 1), respectively. For example, the first rule in category 
Ic of Figure |2] should be interpreted to mean that Si.j,o(^ + 1) — empty 
if Si-ijfi{t) = empty, Si_ij,i(<) = unknown-parity, Si_ij_i,i(i) = 
odd-normal or even, Si-ij-i,o(0 ~ 1' ^^'^ Si,j-i,o(0 = ®- 

A rule in category Id, le, or If is represented as an L-shaped collection 
of three cells, corresponding to the configuration in Figure lb Specifically, 
the upper-right square in each rule corresponds to cell A in Figure lb that 
is, the cell that is evolving, and the color of the square denotes the future 
state into which that cell evolves given the state in its neighborhood. If this 
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cell has coordinates (z, j, 1), the cells to its left, right, and bottom correspond 
to current states of cells G{i + 1) and H{i,j — 1, 1), respectively. For 
example, the first rule in category le of Figure |2] should be interpreted to 
mean that Si,j^i{t + 1) = odd-normal if Sij+i,i(t) = odd-normal and 
Sij,o(0 8 or 2. 

Rule 1. If a cell's neighborhood corresponds to that of any of the evolution 
laws (Figure^, then its state becomes the state of the cell in that evolution 
law corresponding to cell A. 

Rule 2. If the neighborhood of a cell in the top layer does not correspond to 
any evolution law, then that cell's state becomes unknown-parity. 

Rule 3. If the neighborhood of a cell in the bottom layer does not correspond 
to any evolution law, then that cell's state becomes empty. 




FIGURE 3: CAl in action for T" = 7 

In Figure [3] if we treat the non-empty cells in each row as ternary digits, 
the CA gives a sequence of iterates: 

7, 11, 17, 26, 13, 20, 10, 5, 8, 4, 2, 1, 2, 1, 2, 1. 

As we will prove, CAl thus mimics the following modified Collatz map Ti : 
IN ^ IN defined by 

/(3n + l)/2, n=l mod 2 
\n/2, n = mod 2. 

Theorem 1. For every m, there exists some i so that after i iterations of 
CAl, the states of cells not in the state of empty in row k of the bottom 
layer correspond to the ternary digits ofTi{n), where n is the input, for all 
<k <m. 
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Proof. In base three, parity of a number is the same as the parity of the sum 
of that number's digits. So whereas the bottom layer of the CA represents 
transformations of the iterates, the upper layer serves to add digits (mod 2) 
to determine parity. As we shall show, the latter goes about doing this by 
computing successive partial sums, starting with the most significant digits. 

In other words, given an iterate X]fc=o^ i^dk, after i generations have 
elapsed for the CA, the cell above the cell representing djv-j+i will have 
state equal to the parity of '}2k=N-i+i ^k, so after iV + 1 generations, the 
state of the cell above will be the parity of the iterate. 

We claim that the rules in categories Id, le, and If are exhaustive and 
properly determine the parity of an iterate in this way. We will denote cell 
(i, j) to be the cell that is evolving. 

If Si j+i i(i) ~ unknown-parity, the CA is determining the parity of the 
leading digit 0), so Sij.i(< + 1) is simply the parity of Si,jfi{t), which 
can take on three possible values-0, 1, 2-giving the three evolution laws in 
category If. 

Otherwise, if Si.j,o(^) 7^ empty, the CA is calculating the parities of partial 
sums of the k significant digits, and there are three possible states for 0)- 
0, 1, 2- and two for + 1, l)-odd-normal, even-making for the six 
evolution laws in categories le. In each case, Sij^i{t + 1) = Si.j_i.i(t) + 
Si,j.o{t) in the sense of adding parities so that the sum of state even and 1 is 
odd-normal and the sum of odd and 1 is even, for example. 

If Sijfi{t) = empty, the CA has reached the right end of the iterate. If the 
iterate is odd and the coordinates of the units digit are {x, y, 0), the CA sets 
the state of {x, y — 1, 1) to odd-special to signal the CA to append a 1 to 
the iterate (see below), giving the first rule in category Id. The last two rules 
in this category tell the CA to leave all top-layer cells in row i to the right of 
{i, 1) (for odd iterates, to the right of 1)) in state unknown-parity. 

Therefore, the rules in categories Id to Ig properly determine the parity of 
iterates. 

It suffices to show that the rules in categories la to Ic are exhaustive and 
properly mimic the Collatz map in the bottom layer. In particular, we will 
show that if T^{n) — X]i=o '^^^i ^'^^^ and cell (x, y + i, 0) has state di for 
< i < Ar,thenifT'=+i(n) = E^'o ^'^^j, cell {x+l,y+i,Q) will eventually 
have state d'i for < i < N' . lfT^{n) is odd and ifT^+'^{n) = Y.'^^^ 3'd'„ 
cell (cc + 1, y + i - 1, 0) has state d'^ for < i < N'. 

In our proof, we will denote the cell in the bottom-left of each evolution 
law in Figure|2] that is, the cell that is evolving, by (z, j, 0). We will continue 



6 



to denote the current iterate by X^fcLo ^'^'^'fe iterate in the preceding 

row by Y.k=o ^^'^k- 

If Si-i.j-i.o(i) = empty, the CA is determining the rightmost digits of the 
iterate on row i. The first and last rows of category lb in Figure [2]represent 
the two possible cases for (si_i.j i(t), Si-ij-i,i(t)) for when the previous 
iterate is even or odd, respectively. In the former case, the CA should mimic 
division by two so that 2d'i^ = 2si^jfi{t + 1) = Si_i j-.o(0 ~ do (mod 3). In 
the latter case, the CA should mimic r(n) = (3n + l)/2. When Si_i = 
odd-special as in the last rule of category lb, the CA effectively treats the 
empty cell {i — 0) as being in state 1 so that 2dQ = 2sijfi{t + 1) = 1 
(mod 3), meaning the state of {i, j, 0), the units digit of the new iterate will 
be 2. The digit d'l must then satisfy 2{3d'i + cZq) = 3(io + 1 (mod 9), giving 
the second row of rules in category lb. 

If Si-ijfi{t) — empty, the CA is determining the leftmost digits of the 
iterate on row i. By the above discussion of the process by which the CA 
determines parity, cell {i — 1) must be in state unknown-parity. Fur- 
thermore, Sj j_i,o(^) cannot exceed Si_i.j_i_o(t)- Otherwise, if we denote x 
to be the iterate in row i — 1 and y the iterate in row i, then if J2k=o ^'^c^fc is 
even, x/2 < y, and if it is odd, {3x + l)/2 < y, both contradictions! In fact, 
for there to be no such contradictions when Si_ij_io(i) is 1 or 2, Sij-ifi{t) 
must be strictly less than Si_i o(i)- This gives exactly one possibility for 
each choice of Si-i.j-i,o(i) from {0, 1, 2}, and in each case, it is clear that 
because — 1,0) represents the leftmost digit, Si,j^o{t + 1) = empty, 
giving the rules in category Ic. 

If none of the cells in the neighborhood of are empty, the CA is de- 
termining the inner digits of the iterate, so s,i_i j i(i) and Sj_i i(t) must 
be even or odd-normal as in the rules in category la. First note that multi- 
plication by two in base three involves carries of at most one. Therefore, for 
every pair (si_ij^o(i),Si_ij-i,o(t)) G {Q , 1 , 2} x {0 , 1 , 2}, 2sij-ifi{t) = 
■Si-i j-i,o(i) + £ (mod 3), where e E {0, 1}, giving the 18 possible neigh- 
borhoods in category la. In each of these cases, j + 1) satisfies 

2 (3sij,o(< + 1) + Sij-i,o(t)) + e = 3sj_ij^o(i) + Si-i,i-i,o(0 (mod 9) 
giving the evolution laws in category la. □ 

3 CA2-BASE4 

This cellular automaton acts on a two-dimensional infinite grid of square cells 
each in one of five possible states — represented in Figure [5] by the corre- 
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spending colors in parentheses — 8 (white), 1 (light gray), 2 (dark gray), 3 
(black), and empty (white). Additionally, states through 3 each have an 
additional attribute of either -odd or -even, indicated by a circle or a square 
within the cell, respectively (e.g. a state of 0-odd is a white cell containing a 
circle). 

Let the rows be indexed top to bottom by IN U {0} and the columns right 
to left by and denote the cell in the ith row and jth column by {i, j) and 
its state in time period t by Sij {t). 

As shown in Figure|4] the neighborhood of cell (i, j) (white) is the set of 
cells - 1), {i - {i - IJ - 1)} (dark gray). 

FIGURE 4: Neighborhood of CA2 

If input n has quaternary representation X^iLo ^''^^ smallest i for 

which di 7^ is L, then the CA is initialized as follows: for an arbitrary k 
and i from L to A^, we have so.fc+j(0) — d,i ™d Sx,y{0) — empty for all 
(x,2/)^{(0,fc),...,(0,fc + Ar)}. 

In Figure |5] each 2x2 grid represents one evolution law, with each cell 
playing the same role as it does in Figure |4] In other words, the bottom- 
left cell in each 2x2 grid represents the cell that is transforming, and its 
color represents the future state into which it evolves given the states of its 
neighborhood. The colors of the remaining three cells represent the current 
states of the neighborhood that determine the state into which the bottom-left 
cell evolves. For example, the first of eight evolution laws in category 2b 
should be interpreted to mean that Sij {t+1) — empty if Si^ij (t) — empty, 
Si_ij_i(t) — 1-odd and Sij^i{t) = 3-even. 

Rule 1. If a cell 's neighborhood corresponds to one of the evolution laws 
(Figure^, then that cell's state becomes the state of the cell as described in 
the evolution law. 

Rule 2. If a cell's neighborhood does not correspond to any evolution law, 
then that cell 's state becomes empty. 

In Figure [6] if we treat the non-empty cells in each row as quaternary 
digits, CA2 gives a sequence of iterates: 

7,22,11,34,17,13,5,1,1. 
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FIGURE 5: CA2 Evolution laws 



As we will prove, because CA2 mimics the following modified CoUatz map 
T2 : IN ^ IN defined by 




(37i + l)/4'=(3»+i), ^^-^ 2 



where 4'^'^^^ is the largest power of four dividing x, by eliminating trailing 
zeroes. 

Theorem 2. For every m, there exists some i so that after i iterations of 
CA2, the states of cells not in the state of empty in row k correspond to the 
quaternary digits ofT2{n), where n is the input, for all < k < m. 
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FIGURE 6: CA2 in action for = 7 

Proof. Although parity in base four is given by the parity of the units digit, 
each digit must be represented by two states depending on the parity of the 
iterate for the CA to know whether to mimic T{n) = n/2 or Tin) = 3n + 1. 

The former function is equivalent to T{n) = 2n/4 and thus to multiplying 
by two and removing the resulting trailing zero; the latter function is equiva- 
lent to T{n) = {An + 1) — n and thus, roughly speaking, to subtracting the 
iterate from a copy of itself shifted left one place and right-appended with a 
digit of 1 . 

We will denote the cell in the bottom-left of each evolution law in Figure|5] 
that is, the cell that is evolving, by (i, j), the current iterate by J2'k=o ^'^'^Jc 
and the iterate in the preceding row by X^fcLo 4*^6?*;. 

We first prove that the rules in categories 2a to 2c are exhaustive in mim- 
icking T{n) = n/2. Note again that multiplication by two involves carries of 
at most one. 

If Sj j_i(t) = empty as in the rules in category 2c, the CA is determin- 
ing the rightmost digits of the iterate on row i. By design, do is necessarily 
nonzero, so for the iterate in row i — 1 to be even, c?q = 2, giving the first rule 
in category 2c. Here, Sij{t + 1) = 2si-ij{t) (mod 4). There are then four 
choices- 0-, 1-, 2-, 3-even- for di, represented by Si_i ,, (t) in the last four 
rules of category 2c, and because there is a carry of one from multiplying do 
by 2, df) = Sij{t + 1) = 2s,;_ij + 1 (mod 4). 

If Si_i j (t) — empty as in the rules in category 2b, the CA is determining 
the leftmost digits of the iterate on row i. If Si-i,j-i{t) ^ empty, there are 
three possible states for cell {i — 1, j — 1), namely 1-, 2-, 3-even, and in 
each of these cases, there are Si j^i — 2si_i_j_i + e, where e G {0, 1}, 
giving the first six evolution laws in category 2b. In these rules, the state into 
which cell (i, j) evolves is chosen so that Asij + Si.j_i — 2si_ij_i(i) + e. 
If dN' ~ Sij — 1, then the cell to the left of this should be empty to mark the 
end of the iterate, giving the last rule in category 2b. 
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If none of the cells in the neighborhood of are in state empty as in 
the rules in category 2a, the CA is determining the inner digits of the iterate. 
Cells {i — and (i — 1, j — 1) can be in any of four states, namely 0-, 1-, 

2- , or 3-even, and for each of these 16 cases, Si.j-i — 2si_i j_i + e where 
e S {0, 1}, giving the 32 evolution laws in category 2a. In these rules, the 
state into which cell evolves is chosen so that 

4sjj + Sjj_i = 2 (4s,_i,j + Si-ij-i) + 1 (mod 16). 

We next prove that the rules in 2d to 2h are exhaustive in mimicking 

T(n) = 3n + 1. 

If Sij-i{t) = empty as in the rules in category 2f, the CA is determining 
the rightmost digits of the iterate on row i. Because the iterate on row i — 1 
must be odd, there are two choices for do- 1-odd or 3 -odd- represented by 
cell — 1) in the two rules of category 2f's right column (Figure |5|l. If 
do = 3-odd, Sdo + 1 = 2 (mod 4) so that dg = ,,(< + 1) = 2-even, but 
if do — 1-odd, 3do + 1 = (mod 4) so that Sij{t + 1) = empty because 
CA2 deletes trailing zeroes. 

In the former case in which cIq has been determined, the CA then turns to 
evolution laws in category 2d to determine the inner digits of the iterate on 
row i. In the latter case in which cI'q has not been determined, there are four 
possible values for di, namely 0-, 1-, 2-, 3-odd, giving the four evolution 
laws of category 2f's left column. By our interpretation of 3a;+l as (4a:;+l) — 
X, the state into which cell (z, j) evolves is chosen so that cIq ^ Sij{t + 1) = 
Sj_ij_i(t) — Sjj (i) (mod 4) if Si^j (t+1) ^ 0-odd; otherwise, Sij(i+1) = 
empty because CA2 deletes trailing zeroes. Proceeding as above, the CA 
turns to evolution laws in category 2d to determine the iterate's inner digits if 
do has been determined; otherwise, the above process is repeated until dp is 
found. 

If Si-ij{t) — empty as in the rules in categories 2e and 2h, the CA is 
determining the leftmost digits of the iterate on row i. Observe that beyond 
the parity of this iterate, the rules in these categories are the same. 

In this case, if Si_i j_i(t) empty, {i — l,j — 1) can be in any of 
three different states, namely 1-, 2-, 3-odd, and — 1) can be any base- 
four digit in either parity. Certainly the state into which evolves, if 
not empty, must have the same parity as — 1) (e.g. if Sij_i(<) = 

3- even, Sij{t + 1) = n-even for some digit n). In particular, Sij{t + 1) = 
Si_ij_i(t) — Si-ij{t) — e, where e equals 1 if Sij-i + Si_ij_i > 4 and 
otherwise (see discussion of category 2d and 2g rules below). 

On the other hand, if Si_i j_i(i) — empty, (i — 1, j) must be in state 
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empty as well, and can be any nonzero digit in either parity. Here 

Sij-i{t) must be d'j^,: otherwise, if we denote (n) to be the iterate in 
row i — 1 and y to be that in row i, we have that y > T'^+^(n). Therefore, 
Si,j{t + 1) = empty. 

If none of the cells in the neighborhood of (i, j) are in state empty as in 
the rules in category 2d and 2g, the CA is determining the inner digits of the 
iterate on row i. Observe again that beyond the parity of this iterate, the rules 
in these categories are the same. Each cell in the neighborhood can be any 
base-four digit, and the state of — 1) can have either the odd or even 
attribute, giving 64 evolution laws for both category 2d and 2g. Because 
3x + 1 = (4a:: + 1) - x, Sij-i{t) = Si_ij_2(^) - ~ e (mod 4), 

where e e {0, 1}. 

Certainly if Sjj_i(i) + Si-ij-i{t) > 4, the subtraction Si-i,j-2{t) — 
Sj_i j_i — e was performed with a borrowing from Si-i.j (t), so in this case, 
Sij(t + 1) = Si_ij_i(t) - Si_ij(t) - 1 (mod 4). 

If e = 1 and Sij^i{t) + Si_ij_i(f) = 3, then Si_i,j_2(i) = Q- Call 
this digit dk of the iterate on row i — 1. Because e = 1, dk-i = 8 and there 
had to be a borrowing so that dk~2 = ®- Continuing thus, we arrive at a 
contradiction, as the units digit of the iterate on row i — 1 is nonzero. 

So if Sij-i{t)+Si-ij-i{t) < 4, the subtraction Sj_i.j_2(^)— Sj-i j-i — e 
was performed with no borrowing from Si^i,j{t), meaning Sij{t + 1) = 
Sj_ij_i(t) - Sj_ij(t) (mod 4). □ 

4 CA 3 - BASE 2 

This cellular automaton, most efficient and most easily implementable CA, 
acts on an infinite two-dimensional grid of square cells, each of which can 
take on one of three states — ^represented in Figure |8] by the corresponding 
colors in parentheses — 8 (dark gray), 1 (light gray), or empty (white). The 
initial state of the CA is defined to be such that all cells are in the state of 
"empty." 

Again, let the rows be indexed top to bottom by M U {0} and the columns 
right to left by Z, and denote the cell in the ith row and jth column by {i, j) 
and its state in time period t by Sij (t). 

As shown in Figure|7] the neighborhood of cell (i, j) (white) is the set of 
cells - 1), (z - {^ - 1,J - 1), {^ - 1,J - 2)} (dark gi'ay). 

If input n has binary representation X^ilo ^ '^i smallest i for which 

di 7^ is L, then the CA is initialized as follows: for an arbitrary k and i 
from L to N, we have so,/c+i(0) = di and Sx^y{0) = empty for all {x, y) ^ 
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FIGURE 7: Neighborhood of CAS 



{{0,k),...,{0,k + N)}. 



3a 



3b 



3c 



3d 



FIGURE 8: CA3 Evolution laws 

In Figure |8] each group of five cells represents one evolution law, with 
each cell playing the same role as it does in Figure |7] In other words, the 
bottom-left cell in each rule represents the cell that is transforming, and its 
color represents the future state into which it evolves given the states of its 
neighborhood. The colors of the remaining four cells represent the current 
states of the neighborhood that determine the state into which the bottom- 
left cell evolves. For example, the first of six evolution laws in category 3c 
should be interpreted to mean that Si,j{t + 1) = if Si-ij{t) = empty, 
Si_ij_i(t) = 1 and Sij-i(<) = Q, and Si_ij_2(t) = Q- 

Rule 1. // a cell 's neighborhood corresponds to one of the evolution laws 
(Figure^, then that cell's state becomes the state of the cell as described in 
the evolution law. 

Rule 2. If a cell's neighborhood does not correspond to any evolution law, 
then that cell's state becomes empty. 

In Figure |9] if we treat the non-empty cells in each row as quaternary 
digits, CA3 gives a sequence of iterates: 

7,11,17,13,5,1,1. 

As we will prove, because CA3 mimics the following modified Collatz map 
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FIGURE 9: CA3 in action for T° 



Ti : M ^ F defined by 




(37i + l)/2'=(3»+i), j^Q^ 2 

n/2, n = mod 2, 



where 2'^^^^ is the largest power of four dividing x, by ehminating traiUng 
zeroes. 

Theorem 3. For every m, there exists some i so that after i iterations of CAS, 
the states of cells not in the state of empty in row k correspond to the bits of 



Proof. Like in base four, parity in binary is given by the parity of the units 
digit, but unUke CA2, CA3 can mimic T{n) = n/2 simply by representing 
trailing zeroes in an iterate as cells in the empty state. What remains is for 
the CA to mimic T{n) = 3n + 1. 

This function is equivalent to T{n) = {2n + 1) + n and thus, roughly 
speaking, to adding the iterate to a copy of itself shifted left one place and 
right-appended with a digit of 1 . 

We will denote the cell in the bottom-left of each evolution law in Figure|8] 
that is, the cell that is evolving, by the current iterate by X^^o ^'^'^Jc 

and the iterate in the preceding row by J2k=o 2'^c?fe- 

We will prove that the rules in categories 3a to 3d are exhaustive in mim- 
icking T{n) = 3n + 1. Again, multiplication by two involves carries of at 
most one. 

If Sij-i{t) = empty but Si^ij{t) ^ empty as in the rules in category 
3b, the CA is determining the rightmost digits of the iterate on row i. By 
design, do must be 1, giving the first rule in category 2c. Here, Si,j(i, + 1) = 
3sj_i + l (mod 2). There are then two choices- 0, 1- for c?i, represented 
by Si-ij{t) in the next two rules in category 2c's left column. In these two 
cases, the state into which evolves is chosen so that 3{2di + do) + 1 = 
3{2si^i,j{t) + Si^i,j-i{t)) = Sij{t + 1) (mod 2) if Sij ^ and, because 
CA3 treats trailing zeroes as empty cells, empty otherwise. 



3 



(n), where n is the input, for all < k < m. 
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In the former case in which the position of d[-^ = Si,j{t + 1) has been 
determined, the CA then turns to evolution laws in category 3a to determine 
the irmer digits of the iterate on row i. In the latter case in which has not 
been determined, there are two possible values for di, namely 8, 1, giving the 
bottom two evolution laws in category 3b. By our interpretation of 3a; + 1 as 
{2x + 1) + a;, the state into which cell {i, j) evolves is chosen so that 

d'o = Si^j{t + 1) = Si_ij_i(i) + Sj_i,j(i) + e (mod 4) 

where e = 1 if Si-\j-\{t) + Si-\^j-2{t) > Sij-i and e = otherwise (see 
discussion of category 3a rules below) , provided that Sij{t + 1) 7^ 0; if not. 
Si J {t+l) = empty because CA3 treats trailing zeroes as empty. 

Again, in the former case in which the position of d'^ has been determined, 
the CA then turns to rules in category 3a to determine inner digits. In the 
latter case in which d'^ has not been determined, there are two possible values 
for d2, namely 0, 1, giving the top two evolution laws in category 3b's right 
coluinn. The state into which {i, j) evolves, either empty or 1, is determined 
in the same way as above. One more, the CA turns to evolution laws in 
category 2d to determine the iterate's inner digits if d'^ has been determined; 
otherwise, the above process is repeated until the position of dg is found. 

If Si-ij{t) = empty as in the rules in categories 3c and 3d, the CA is 
determining the leftmost digits of the iterate on row i. 

In this case, if Sj_ij_i(t) ^ empty, it represents dN and thus must be 
1. There are only two cases in which at least one of cells (i — 1, j — 2) and 

— 1) is empty: i) the iterate on row i — 1 is of the form x = J2h=o 
so that T{x) = 4*^+^ and the iterate on row i is thus 1, ii) the iterate on row 
i — 1 is a; = 1 so that the iterate on row i is also 1. These cases give rise to the 
first two rules of category 3d, respectively, where Sij{t + 1) is empty rather 
than Q because CA3 deletes trailing zeroes. The last rule of 3d then accounts 
for the case in which the iterate on row i is 1 . 

Beyond these two cases, if neither cell {i — 1, j — 2) nor cell {i,j — 1) is 
in the empty state, there are two possible states for each, giving the first four 
rules in category 3c. Here, the state into which {i, j) evolves is again chosen 
so that 

Sij{t + 1) = Si_ij_i(t) + Si-i,j(t) + e = Sj_ij_i(i) + e (mod 4) 

where e = 1 if Si_ij_i(t) + Si-ij-2{t) > Si,j-i and e = otherwise. 

On the other hand, if Si-ij-i{t) = empty, Si_i j_2(i) represents dN- 
There are two possible non-empty states for cell {i, j — 2)- 8, 1- giving the 
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last two rules in category 3c. In the former case, must represent rf'^,; 
otherwise, if we denote this iterate to be y and the preceding iterate to be x, 
y > T3{x), a contradiction. Therefore Sij{t + 1) = 1. In the latter case, 
— 1) must represent d'pf, for the same reason, so Sij + 1) = Q. 
Finally, if none of the cells in the neighborhood of are in state empty 
as in the rules in category 3a, the CA is determining the inner digits of the 
iterate. Each of the cells in the neighborhood can be in any of two possible 
states, giving the sixteen rules in this category. Because 3x+l = {2x+l) — x, 
Sij-i{t) = Sj_ij_2(t) + + e (mod 2), where e e {0, 1}. 

Sij-i{t) < Sj_i,j_2(t) + Si_ij_i(i) < Sj_ij_2(t) + Si_ij_i(i) + e 

impUes the addition Si^ij-2{t) + Si-ij-i{t) was performed with a carry 
over to so Si,j{t + 1) = Si-i,j-i(t) + Si_i.j(t) + 1 (mod 2). On 

the other hand, Sij-i{t) > Si-ij-2{t) + Si-ij-i{t) implies there was no 
carry, and Sij{t + 1) = Si_ij_i(t) + Si-ij{t) (mod 2). □ 

5 CONCLUSION 

Call the n-efficiency of any of the above CAs the number of iterates that have 

been computed by the CA upon reaching 1 given input n, divided by the total 
stopping time of that input. The average n-efficiency over all n such that 
1 < n < 2^^ was determined to be roughly 69.4%, 63.7%, and 32.2% for 
CAl, CA2, and CA3, respectively. 

In addition, observe that it is not necessary to initialize any of these CAs 
with only one input. In fact, these CAs are able to compute the trajectories of 
an arbitrary number of inputs at a time. As an example, in the case of CA3, 
let the set of inputs be {nj}i<j<M, where input rij has binary representation 
X^^o ^^'^i ^^^^ that the smallest i for which 7^ is some Lj. The CA is 
initialized as follows: pick an arbitrary fco € ^- For all j such that 1 < j < 
M, let kj = kj_i + Nj_i + ej, where ej is chosen so that the trajectories 
of 'J2i'=o "^^^i 12i=^a^ do not collide on the grid. Then for each j 

such that 1 < j < M, letso,fc,+i(0) = for < « < Nj. Let Sx,y{0) = 
empty for all other pairs (x, y). 

Note however, that the difficulty in initializing a CA to compute multiple 
trajectories in parallel this way lies in the unpredictability of the behavior of 
ej: it is not certain how large the space between two inputs must be for their 
trajectories not to colUde. 

As a workaround, if we place each iterate whose trajectory we want to 
compute on its own grid and, roughly speaking, stack those grids on top of 
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FIGURE 10: Parallel computing for initial iterates n =183, 120767, 
53132499 

one another to form an arbitrarily high three-dimensional grid of cubical cells, 
then this CA can verify multiple inputs in parallel. 




FIGURE 1 1 : Parallel computing in three dimensions 

In conclusion, we have found three CAs that mimic the behavior of the 
CoUatz map in bases two, three, and four, and streamline the process of veri- 
fication of the CoUatz Conjecture by bypassing calculation of certain parts 
of the trajectory. Beyond speeding the progress towards potentially find- 
ing a counterexample, this work affords insights into the distribution of the 
(mod 2*^) residues of iterates, which may justify the existing heuristic prob- 
abilistic argument that trajectories tend to decrease [7]. 

Topics for further study include applying the parallel computing model we 
propose to similar iterative computations, exploring whether it is possible to 
create evolution laws that depend entirely on digits from the previous iterate, 
and exploring other potentially viable radices for new CAs. With regards 
to the latter, it is known that the distribution of (3/2)'^ (mod 1) is closely 
related to the "sorting properties" of the CoUatz map, that is, whether the 
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function distributes odd iterates equally among the residues (mod 2*^) [3]. If 
no reasonable CA exists in base 3/2, it could still be worth exploring possible 
CAs for computing [(3/2) J . Implementations in Java of CA2 and CA3 can 
be found at www . sitanchen . com/collatz4 . htrnl] and |www.sitanchen.| 



com/ collatz . html respectively. 
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